First tests

import imageio.v3 as iio
import numpy as np

hm = iio.imread("Obsidian.png", mode="F")
hm = (hm - hm.min()) / (hm.max()-hm.min()) #normalize heightmap

hm

array([[0.46236876, 0.4609779 , 0.4617027 , ..., 0.45954788, 0.46223164,
        0.46221203],
       [0.4605861 , 0.46070364, 0.46031186, ..., 0.4578436 , 0.45868593,
        0.4593128 ],
       [0.45790237, 0.45735386, 0.4593128 , ..., 0.4560022 , 0.45794153,
        0.45737344],
       ...,
       [0.4639163 , 0.4648566 , 0.46609074, ..., 0.46571854, 0.4631915 ,
        0.46417096],
       [0.46325028, 0.46407303, 0.46491536, ..., 0.46407303, 0.46350494,
        0.46197698],
       [0.46236876, 0.46266258, 0.4624471 , ..., 0.46250588, 0.462976  ,
        0.46217287]], dtype=float32)

import trimesh

mesh = trimesh.load("test.stl")

mesh.show()